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1.0  INTRODUCTION 
1.1  PROGRAM  BACKGROUND 

Extensive  work  has  been  devoted  to  the  development  of  infrared  detection  systems  for  a 
variety  of  orbital  surveillance  functions.  Various  thrusters  are  used  aboard  satellites  for 
attitude  control  and  station-keeping  purposes,  and  interest  in  condensation  phenomena 
within  the  plumes  of  these  thrusters  was  stimulated  by  the  realization  that  gaseous  plume 
constituents  could  condense  into  conglomerates  of  larger  particles  which  could  ultimately 
pass  within  the  field  of  view  of  the  infrared  detection  system.  These  conglomerates  could 
affect  the  satellite-borne  sensors  by  several  mechanisms: 

1.  The  conglomerates  can  a)  emit  or  b)  absorb  electromagnetic  radiation  within  the 
spectral  region  to  which  the  infrared  detection  system  is  sensitive. 

2.  The  conglomerates  can  provide  large  electromagnetic  radiation  scattering 
centers  which  can  scatter  radiation  to  the  detector  from  a  source  normally  out  of 
the  field  of  view. 

As  a  result,  the  sensor’s  background  noise  can  be  increased,  false  target  information  can  be 
provided,  and  the  sensitivity  of  the  detector  can  be  decreased. 

In  FY73  initial  experimental  studies  of  plume  condensation  were  begun  under  the  Air 
Force  Rocket  Propulsion  Laboratory  (AFRPL)  CONSCAT  Program.  The  objectives  of  this 
program  were  to  determine  the  amount  of  Rayleigh  scattered  radiation  from  the 
homogeneous  nucleation  of  nitrogen  (NJ  in  a  supersonic  flow  from  a  conical  nozzle  for  a 
known  incident  radiant  energy  at  a  particular  visible  wavelength.  Infrared  scattering  was  to 
be  inferred  by  using  the  Rayleigh  scattering  relations  (Refs.  1  through  3). 

The  axial  profiles  of  Rayleigh  scattered  intensity  were  in  good  agreement  with  the 
intensity  levels  expected  for  an  isentropic,  uncondensed  N2  expansion  for  low  reservoir 
pressure  and  the  dimensionless  axial  position  x  =  x/D,  where  D  is  a  characteristic  throat 
diameter  for  the  gas  source.  The  scattered  intensity  showed  a  subsequent  sharp  increase  with 
x  at  low-to-intermediate  reservoir  pressures,  which  was  interpreted  as  indicating  that 
condensation  onset  was  occuring.  With  continued  increasing  x,  the  scattering  signal  went 
through  a  maximum  and  then  decayed.  This  behavior  was  interpreted  as  indicating 
condensate  growth  followed  by  cessation  of  the  condensation  process.  An  analytic  method 
was  developed  to  calculate  condensation  and  cluster  growth.  Knowledge  of  the  experimental 
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scattering  results  enabled  selection  of  initial  nucleation  conditions  for  the  calculation  such 
that  the  calculated  properties  of  condensed  Ni  nozzle  flow  such  as  gas  density,  temperature, 
spatial  locations  of  onset  and  cessation  of  condensation,  and  Rayleigh  scattered  intensity 
were  in  satisfactory  agreement  with  the  experimental  results.  Condensate  mass  fractions  of 
10'2  or  less  produced  experimentally  observable  increases  in  the  rotational  temperature  of 
the  gas.  However,  no  noticeable  change  in  the  monomer  number  density  was  detected. 
Estimates  of  condensate  cluster  diameter  ( =  100  A)  supported  the  validity  of  the  application 
of  Rayleigh  scattering  for  the  study  of  the  condensation  process  in  expansion  flows. 

1.2  OBJECTIVES 

It  was  concluded  in  Ref.  1  that  in  order  to  develop  an  a  priori  condensation  calculation 
procedure  for  a  wide  range  of  flow  conditions  one  must  have  knowledge  of  the  scaling  laws 
of  condensation  for  flow  source  geometry,  reservoir  conditions,  and  molecular  parameters. 
Consequently,  the  CONSET  Program  was  formulated  with  the  objective  of  determining 
experimentally  the  onset  and  growth  properties  of  condensate  clusters  in  typical  exhaust 
plume  flow  fields  and  the  dependence  of  the  condensation  process  on  nozzle  geometry, 
reservoir  conditions,  molecular  parameters,  and  flow  composition.  The  program  was 
initially  divided  into  three  task  areas.  The  first  task  was  to  investigate  the  effects  of  variation 
of  reservoir  pressure  and  temperature,  throat  diameter,  and  nozzle  expansion  angle  on  the 
spatial  laws  of  condensation  onset  and  growth  for  Ni  flowfields.  The  second  task  was  to 
extend  the  condensation  scaling  law  study  to  other  species  such  as  argon  (Ar),  oxygen  (Oi), 
carbon  monoxide  (CO),  nitric  oxide  (NO),  carbon  dioxide  (CO2),  ammonia  (NH3),  water 
vapor  (H2O),  and  hydrogen  chloride  (HC1),  as  well  as  some  binary  mixtures  of  these  species. 
The  third  task  was  to  investigate  the  condensation  process  in  both  an  actual  thruster  plume 
and  a  simulated  thruster  plume.  The  actual  thruster  to  be  used  was  a  0.1-lbf 
monopropellant  hydrazine  engine. 

The  experimental  investigations  were  conducted  using  four  noninterfering, 
nonperturbing  flow  diagnostic  techniques.  Laser-Rayleigh  scattering  was  used  to 
characterize  the  distribution  and  growth  of  clusters  as  well  as  condensation  onset.  Laser- 
Raman  scattering  was  used  for  measurement  of  plume  species  number  density  and 
temperature  in  the  higher  density,  near-field  region  of  the  expansion,  and  electron  beam 
fluorescence  was  used  for  measurement  of  N2  gas  density  and  temperature  in  the  far  field  of 
the  expansion.  The  mass  spectrometer  probe  was  used  for  far-field  measurements  of  the 
monopropellant  thruster  relative  species  concentrations. 
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1.3  RAYLEIGH  SCATTERING  TECHNIQUE 

The  basic  equations  and  rationale  for  the  application  of  Rayleigh  scattering  to  the  study 
of  condensing  gas  flow  fields  have  been  given  in  detail  in  Refs.  4  and  5.  Therefore,  only  a 
brief  summary  need  be  given  in  this  report. 

For  an  incident  laser  beam  of  wavelength  \q  and  intensity  I0  focused  within  a  flowing 
pure  gas  sample  of  number  density  n  with  species  polarizability  a,  the  scattered  intensity,  I, 
which  is  normalized  by  I0,  is  given  by 

I  =  Kn  a2/**  (1) 


in  which  K  is  a  coefficient  containing  transmission  and  calibration  factors.  For  a  scatterer  of 
radius  “a”  which  is  characterized  by  bulk  properties,  it  is  known  that  a  is  proportional  to 
a3,  which  indicates  the  sensitivity  of  Rayleigh  scattering  to  scatterer  size.  Assuming  the 
condensing  flow  field  to  be  composed  of  a  collection  of  gas  phase  monomers  and  molecular 
clusters,  or  i-mers,  where  i  is  the  number  of  molecules  per  cluster,  the  single  Rayleigh 
scattering  intensity  with  polarization  vector  parallel  to  the  incident  beam’s  plane  of 
polarization  is 


I'd)  = 


(2) 


where  nQ  is  the  reservoir  number  density  of  the  flow  field.  The  scattered  intensity  I' (II) 
includes  the  further  normalization  provided  by  the  scattered  intensity  from  a  collection  of 
monomers  of  number  density  nD- 


For  an  uncondensed,  isentropic  expansion. 


!'(»)  = 


1°  ( n ) 


(3) 


Super-  and  subscript  zeros  denote  isentropic  and  reservoir  conditions,  respectively.  The  axial 
variation  of  I°(ll)  is  provided  by  the  method-of-characteristics  solution  (MOCS)  (Ref.  6)  for 
nozzle  flow  and  by  the  Sherman-Ashkenas  theory  (Ref.  7)  for  sonic  orifice  flow. 
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Deviation  of  the  measured  I  '(H)  from  I0(B)  indicates,  for  these  studies,  the  existence  of 
condensation.  A  direct  measure  of  the  existence  of  clusters  within  the  flow  is  given  by  the 
scattering  function,  f,  which  is  written  as 

f,/_LU! i\_]  (4) 

\r*  ni )/ 

Although  the  scattering  function,  f,  is  an  ambiguous  measure  of  the  simultaneous  increase  in 
the  mean  cluster  size  and  condensate  mole  fraction,  the  axial  variation  of  f  as  a  function  of 
reservoir  conditions  and  nozzle  geometry  yields  empirical  condensation  scaling  laws. 

1.4  RAMAN  SCATTERING  AND  ELECTRON  BEAM  FLUORESCENCE 

Laser-Raman  scattering  and  electron  beam  fluorescence  were  the  diagnostic  techniques 
used  for  plume  static  property  measurements  in  the  near  and  far  fields,  respectively.  The  use 
of  the  electron  beam-induced  fluorescence  radiation  resulting  from  inelastic  electron- 
molecule  collisions  for  the  measurement  of  N2  gas  density  and  temperature  in  flow  fields  has 
been  discussed  and  demonstrated  in  detail  in  Refs.  5  and  8.  Briefly,  a  high-energy  beam 
of  electrons  produces  excited  electronic  states  upon  impact  with  N2  molecular  species.  By 
radiative  decay  these  states  produce  a  fluorescence,  and,  assuming  no  collisional  quenching 
effects,  the  intensity  of  the  fluorescence  is  directly  proportional  to  the  N2  gas  density.  The 
electron  beam  technique  is  normally  used  at  total  gas  number  densities  of  less  than  1015cm  3 
(Ref.  9)  to  avoid  collisional  quenching  effects,  and  for  that  reason  the  technique  is  general¬ 
ly  a  far-field  plume  diagnostic.  The  molecular  fluorescence  consists  of  vibrational-rotational 
band  structure  within  the  electronic  transition  systems,  which,  upon  spectral  dispersion, 
yields  the  temperature  of  the  rotational  energy  mode  (Ref.  8). 

As  a  result  of  nonresonant,  inelastic  collisions  between  laser  beam  photons  and  gas 
molecules,  radiation  is  scattered  from  the  beam  at  frequencies  other  than  that  of  the  laser 
light.  This  phenomenon,  spontaneous  Raman  scattering,  is  discussed  in  detail  in  Refs.  10 
through  13,  and  details  of  application  to  measurement  of  species  number  density  and 
temperature  are  fully  discussed  in  Ref.  14.  Briefly,  a  high-energy  laser  beam  interacts  with 
gas  molecules  and  causes  changes  in  the  rotational  and  vibrational  energy  levels.  As  a  result, 
the  Raman  scattered  radiation  consists  of  molecularly  specific  pure  rotational  lines  as  well  as 
vibrational-rotational  band  structure  from  which  rotational  and  vibrational  temperatures 
can  be  determined.  The  intensity  of  the  lines  and  bands  is  also  directly  proportional  to  the 
molecular  species  number  density.  Because  it  is  a  scattering  process,  the  Raman  intensity  is 
not  subject  to  collisional  quenching  effects;  therefore,  the  Raman  technique  is  useful  in  the 
near-field  plume.  However,  the  low  Raman  scattering  cross  section  ( =  10'28cm2/sr)  generally 
prevents  use  of  the  technique  for  far-field  measurements. 
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1.5  MASS  SPECTROMETRIC  SAMPLING 

The  direct-sampling  mass  spectrometer  probe  developed  at  AEDC  (Ref.  15)  was  used  for 
species  mole  fraction  measurements  on  the  axial  centerline  of  the  monopropellant  hydrazine 
thruster.  In  reality,  the  probe  is  merely  a  miniature  molecular  beam  system  which  uses  a 
quadrupole  mass  spectrometer  as  the  detector.  Immersing  the  probe  in  the  exhaust  plume  of 
a  rocket  engine  forms  a  molecular  beam  of  exhaust  products  which  is  directed  into  the 
sampling  volume  of  the  mass  spectrometer  ion  source.  Subsequently,  ions  are  created, 
extracted,  and  directed  through  a  mass  discriminator  and  then  into  an  electron  multiplier. 
The  result  is  a  set  of  measurable  signals  that  are  proportional  to  the  number  densities  of 
molecules  in  the  exhaust  plume  which  (after  ionization)  possess  particular  charge-to-mass 
ratios.  The  interpretation  of  the  signals  is  dependent  upon  a  reliable  calibration,  a 
knowledge  of  the  additional  daughter  mass  peaks  (cracking  patterns)  created  in  the  ion 
source  from  plume  parent  constituents,  and  the  ability  to  maintain  the  very  stringent 
vacuum  conditions  necessary  to  operate  a  molecular  beam  system. 

2.0  DESCRIPTION  OF  EXPERIMENTAL  SETUP  FOR  SCALING 
LAW  AND  SIMULATED  THRUSTER  STUDIES 

2.1  RESEARCH  VACUUM  CHAMBER 

The  experiments  reported  herein  were  conducted  using  the  4-  by  10-ft  Research  Vacuum 
Chamber  (RVC),  a  stainless  steel  vacuum  chamber  nominally  4  ft  ( 1 .3  m)  in  diameter  and  10 
ft  (3.3  m)  long.  The  chamber,  shown  in  Fig.  1,  is  constructed  in  two  sections:  a  movable 
section  approximately  4  ft  (1.3  m)  long  and  a  stationary  section  approximately  6  ft  (2  m) 
long.  An  additional  spool  piece  1.5  ft  (0.49  m)  long  was  used  to  provide  eight 
instrumentation  ports. 

Initial  pumping  capabilities  were  supplied  by  a  300-cfm  mechanical  pump  for  rough 
pumping  and  a  6-in.-diam  oil  diffusion  pump  with  a  baffle  refrigerated  by  liquid  nitrogen 
for  intermediate  pumping.  Pumping  during  experiments  was  provided  by  a  finned,  gaseous 
helium  cyroliner  at  20  K  with  a  liquid  nitrogen  cryoliner  radiation  shield  located  in  the 
stationary  section  of  the  RVC.  The  blank-off  chamber  pressure  achieved  was  approximately 
10'7  torr. 

2.2  GAS  FLOW  SOURCES 

The  flow  generators  were  sonic  orifices  and  conical  nozzles  attached  to  a  GTE-Sylvania 
heated  source.  This  assembly  was  mounted  on  a  motor-driven,  three-dimensional  traversing 
mechanism  located  in  the  movable  section  of  the  RVC.  The  sonic  orifices  were  1.325,  3.2, 


11 


AEDC -TR-80-16 


and  3.05  mm  in  diameter  with  a  diameter-to-thickness  ratio  greater  than  20.  The  conical 
nozzles  had  a  nominal  throat  diameter  and  nozzle  length  of  1.0  and  5.334  mm,  respectively. 
The  four  nozzle  half-angles  were  14.5,  10.5,  9.0,  and  5.63  deg.  Schematics  of  the  sonic 
orifice  and  conical  nozzle  are  shown  in  Fig.  2.  The  conical  nozzle  for  the  simulated  thruster 
studies  was  fabricated  to  duplicate  the  nozzle  of  the  monopropellant  thruster.  The  half¬ 
angle  was  15  deg,  and  the  throat  diameter  was  0.03  in.  (0.76  mm).  The  length  from  the 
throat  to  the  exit  plane  was  0.356  in.  (9.04  mm),  and  the  expansion  ratio  was  55.3. 

The  gas  reservoir  was  instrumented  with  standard,  calibrated  pressure  and  temperature 
gages.  Gases  were  supplied  from  high-pressure  bottles,  and  two  25.0-nanometer  (nm)  filters 
were  installed  in  the  inlet  line  to  minimize  effects  of  particulate  matter. 

2.3  OPTICAL  SYSTEM 

An  argon-ion  laser  operating  at  a  wavelength  of  514.5  nm  provided  excitation  for  the 
Rayleigh/Raman  scattering  measurements.  As  shown  in  Fig.  1,  the  incident  laser  beam 
polarization  was  rotated  along  the  x-direction,  expanded,  and  focused  onto  the  chamber 
centerline.  Light  scattered  from  the  focal  volume  was  collected  by  an  f/2  lens  system, 
collimated,  and  focused  onto  the  input  slit  of  a  0.85-m,  double-grating  spectrometer.  For 
the  Rayleigh  scattering  measurements,  HN-22  Polaroid1’’  material  was  placed  in  the 
collimated  light  path,  and  a  polarization  scrambler  was  placed  immediately  in  front  of  the 
spectrometer  entrance  slit.  During  scaling  law  studies  the  entrance  slit  aperture  setting, 
collection  optics  magnification,  and  laser  beam  focusing  together  resulted  in  observation  of 
a  1 .5-mm-long,  50-  to  100-gm-diam  cylindrical  scattering  volume.  For  the  simulated  thruster 
measurements  a  2.75-mm  long,  600-gm-diam  cylindrical  scattering  volume  was  observed 
using  an  f/4  collection  lens  system. 

The  radiation  detector  was  a  thermoelectrically  cooled  photomultiplier,  and  the  output 
was  processed  by  an  Ortec-  photon-counting  system  for  either  digital  display  or  strip  chart 
recording. 

For  the  electron  beam  measurements,  the  laser  dump  and  laser  input  aperture  assembly 
shown  in  Fig.  1  were  removed  and  replaced  by  an  RCA  Model  VC2126V4  electron  gun 
assembly  and  cylindrical  Faraday  cup  assembly,  respectively.  The  beam  was  injected  into  the 
chamber  and  flow  field  through  a  1.0-mm-diam  orifice  which  provided  the  necessary 
pressure  drop  to  maintain  the  electron  gun  pressure  at  less  than  6  x  10'6  torr  under  nominal 
chamber  vacuum  conditions.  For  number  density  and  temperature  measurements  the  beam 
current  was  1.0  mA  at  30keV  energy.  The  collection  optics/spectrometer/data  acquisition 
system  was  the  same  as  that  used  for  the  laser  scattering  measurements. 
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3.0  DESCRIPTION  OF  EXPERIMENTAL  SETUP  FOR 
MONOPROPELLANT  THRUSTER  STUDIES 

3.1  RESEARCH  VACUUM  CHAMBER 

Figure  3  is  a  schematic  of  the  monopropellant  thruster  installation  in  the  RVC.  For  the 
thruster  experiments  the  RVC  was  reconfigured  by  using  the  4-ft-long  movable  section  to 
house  a  37-liter  capacity,  liquid  helium-filled  cryopump.  The  section  was  placed  directly  in 
line  with  the  fixed  section  housing  the  cryoliner.  Another  movable  section,  approximately  3 
ft  (0.99  m)  long,  was  connected  to  the  spoolpiece,  and  it  was  in  this  section  that  the  thruster 
was  mounted  on  the  same  traversing  table  as  was  used  for  the  scaling  law  experiments. 

3.2  THRUSTER 

The  thruster  used  in  the  experimental  program  was  a  Hamilton  Standard  REA/CTS 
10-18  monopropellant  hydrazine  thruster  which  employed  a  30-35  mesh  Shell  405  ABSG 
spontaneous  catalyst  with  a  preloaded,  packed  bed  design.  The  exhaust  was  provided  by  a 
conical  nozzle  with  a  0.76-mm-diam  throat  and  an  exit  area  ratio  of  55.  The  thruster  was 
operated  over  a  thrust  range  from  0.44  to  1.10  N  (0.1  to  0.25  lbf)  with  a  nominal  0.14-sec- 
on/9. 86-sec-off  duty  cycle  using  initial  catalyst  bed  temperatures  of  367  K  (200°F),  478  K 
(400°F),  and  589  K  (600°F).  The  engine  was  always  operated  without  its  heat  shield.  A 
complete  description  of  the  engine,  its  properties,  propellant  system,  and  method  of 
operation  may  be  found  in  Ref.  16. 

3.3  DIAGNOSTIC  SYSTEMS 

Exhaust  plume  diagnostic  systems  employed  included  a  mass  spectrometer  probe,  a 
quartz  crystal  microbalance  (QCM),  a  laser  Raman/ Rayleigh  scattering  system,  an  electron 
beam  fluorescence  system,  and  a  particle  collection  network.  The  QCM  and  particle 
collection  network  were  for  providing  information  on  the  contamination  properties  of  the 
exhaust  plume  rather  than  for  condensation  diagnostics;  therefore,  they  are  not  discussed  in 
this  report. 

Figure  4  is  a  schematic  diagram  of  the  special  diagnostic  instrumentation.  Laser 
scattering  excitation  was  provided  by  a  conventional  mode  ruby  laser  with  a  pulse  width  of 
approximately  1  msec.  The  electron  beam  system  was  the  same  as  that  used  for  the  scaling 
law  studies.  A  0.5-m  double  grating  spectrometer  was  used  for  spectral  dispersion  for  both 
laser  scattering  and  electron  beam  measurements.  The  quadrupole  mass 
spectrometer/skimmer  system  was  located  on  the  centerline  at  the  axial  position  450  <  x 
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<650.  The  probe  face  was  conical  with  a  cone  half-angle  of  30  deg,  and  a  15-deg  conical 
half-angle,  pure  nickel  skimmer  with  a  0.020-in. -di am  orifice  was  soldered  to  the  probe 
body.  The  forward  section  of  the  probe  was  cooled  with  gaseous  helium  at  =20  K. 

Mass  spectrometer  data  were  displayed  on  an  oscillograph  as  well  as  acquired  digitally. 
Electron  beam  and  laser  scattering  data  were  acquired  using  a  cooled  photomultiplier  tube 
and  photon  counting.  Data  processing  and  initial  analysis  of  all  data  were  performed  using 
an  online  PDP-8  computer  system  with  line  printer  output. 

4.0  DISCUSSION  OF  SCALING  LAW  DATA 

4.1  RAYLEIGH  SCATTERING  RESULTS  FROM  THE  SCALING 
LAW  EXPERIMENTS 

Figures  5  through  8  show  the  axial  profiles  of  I '(B)  for  the  N2,  O2,  CO,  and  Ar 
expansions  from  sonic  orifices.  Figures  9  through  24  show  the  axial  profiles  of  I' (II)  for  the 
N2,  O2,  CO,  NO,  HC1,  H20,  CO2,  and  Ar  expansions  from  conical  nozzles.  Theoretical 
predictions  as  obtained  from  the  Sherman-Ashkenas  theory  or  the  MOCS  are  also  shown.  It 
is  observed  that  the  onset  of  condensation  is  manifested  by  a  dramatic  increase  of  I '(B) 
relative  to  the  isentropic  prediction.  The  onset  of  condensation  moves  nearer  the  saturation 
point  as  Pu  increases  or  Tn  decreases,  and,  with  the  exception  of  the  region  of  discontinuity 
in  the  MOCS  nozzle  calculations,  I' (II)  is  in  good  agreement  with  the  calculated  values  prior 
to  condensation  onset.  Furthermore,  with  the  exception  of  H20  and  NO,  it  can  be  seen  that 
for  the  lowest  P0  values  or  highest  T0  values  the  metastable  gas  sample  supports  a 
supersaturated  state  for  approximately  30  nozzle  throat  or  sonic  orifice  diameters  before 
condensing.  It  is  also  noted  (generally)  that  the  massive  condensate  growth  for  the  nozzle 
flows  is  rather  abrupt,  whereas  the  massive  condensate  growth  region  for  the  sonic  orifice 
flow  is  preceded  by  a  gradual  deviation  from  the  isentropic  prediction.  The  magnitude  of  the 
Rayleigh  scattering  intensity  in  the  condensate  growth  regions  is  seen  to  increase  rapidly  with 
increasing  P0  and  to  decrease  extremely  rapidly  with  increasing  T0. 

Radial  profiles  of  1' (II)  for  Ni  are  shown  in  Figs.  25  through  27.  In  Figs.  25  and  26  the 
x  =  17.35  and  17.9  positions,  respectively,  are  interesting  in  that  two  scattering  peaks  are 
symmetrically  located  off  the  axial  centerline.  Similar  observations  have  been  reported  by 
Beylich  (Ref.  17)  in  a  study  of  CO2  condensation  in  a  nozzle  flow.  From  Figs.  1 1  and  19,  it  is 
seen  that  condensation  on  the  axial  centerline  has  not  begun  at  x  =  17.35  and  17.9, 
respectively.  However,  Figs.  25  and  26  show  that  onset  has  already  begun  for  r/D,  >  0. 
Photographic  observations  of  the  6  =  14.5-deg  nozzle  flow  field  are  presented  in  Ref.  18, 

and  the  nozzle,  a  dark,  isentropic  expansion  zone,  and  a  bright,  hemispherical  onset  zone 
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are  clearly  evident,  as  is  Alimentary  structure  within  the  condensation  growth  region.  These 
photographs  substantiate  the  radial  profile  results  of  Figs.  25  and  26.  The  radial  profiles  of 
Fig.  27  are  interesting  because  they  demonstrate  the  sensitivity  of  the  condensation  process 
to  reservoir  temperature.  An  increase  of  67  K  in  T0  has  completely  suppressed  the  off-axis 
condensation  at  x  =  17.9.  At  x  =  55.4,  off-axis  peaks  can  again  be  observed,  however. 

Figures  28  and  29  are  axial  variations  of  I'(H)  for  95-percent  N2/5-percent  CO2  and 
90-percent  N2/10-percent  CO2  mixtures,  respectively.  As  demonstrated  in  Ref.  19,  a  simple 
computation  of  reservoir  partial  pressures  and  reference  to  the  pure  gas  axial  variations  of 
I'(l)  for  the  appropriate  nozzle  reveal  that  I ' (H)  has  been  greatly  enhanced  relative  to  the 
level  expected  based  on  the  pure  gas  expansions.  Saturation  locations,  xs,  are  indicated  in 
Figs.  28  and  29  for  both  CO2  and  N2.  Obviously  CO2  is  the  initiator  of  the  onset,  because  N2 
is  not  supersaturated  in  the  region  of  onset.  These  results  demonstrate  the  important  role  of 
mixed  clusters  in  gas  mixture  condensation  processes. 

Figures  30  through  32  demonstrate  the  axial  variation  of  the  scattering  function,  f,  for 
N2  gas  for  conical  nozzles  (Figs.  30  and  32)  and  a  sonic  orifice  (Fig.  31).  These  plots  are 
rather  typical  of  all  the  f-versus-x  plots  obtained  in  the  scaling  law  study  for  those  gases 
which  exhibit  condensation  onset  external  to  the  nozzle.  It  is  noted  that  in  these 
semilogarithmic  plots  the  f  values  for  a  given  P0  (or  T0)  and  nozzle  (or  orifice)  form  straight 
lines,  and  it  is  the  intersection  of  these  straight  lines  with  the  x-axis  that  is  used  to  determine 
the  onset  of  condensate  growth.  These  axial  onset  locations  are  denoted  by  xB.  The  rapid 
increase  in  f  following  onset  is  obvious,  as  are  the  orders  of  magnitude  increase  in  f  as  the 
reservoir  pressure  increases  or  reservoir  temperature  decreases. 

The  vapor  pressure  data  compiled  in  either  Ref.  20  (Hilsenrath  et  al.)  or  Ref.  21  and  the 
isentropic  solution  for  each  particular  flow  field  investigated  are  used  to  obtain  the 
saturation  values  of  pressure  and  temperature,  Ps  and  Ts,  respectively,  as  well  as  similar 
values  at  consensation  onset,  denoted  by  Pg  and  Tg.  The  isentropic  super  saturation  pressure 
ratio,  (sg)°,  is  defined  as 


(’9)’  -  P./P8  <5> 

and  the  isentropic  degrees  of  supercooling,  (Sg)°,  are  defined  as 

W  -  T.-Te  (6) 
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It  should  be  noted  that  the  supersaturation  ratio  defined  here  is  not  the  normal  definition  of 
the  ratio  of  the  pressure  at  onset  to  the  equilibrium  vapor  pressure  at  onset,  Pg/P^g.  The 
definition  of  Eq.  (5)  is  used  because  of  the  unreliability  of  Pvg  values  at  the  low  temperatures 
at  condensation  onset  in  the  expansion  flows  of  these  experiments.  The  (sg)°  values  are 
orders  of  magnitude  lower  than  the  normally  defined  supersaturation  ratios  for  the  flow 
fields  investigated  here.  These  supersaturation  parameters  are  illustrated  in  Fig.  33,  which  is 
a  diagram  of  the  expansion  process. 

Table  1  is  a  tabulation  of  the  saturation,  condensation  onset,  and  supersaturation 
parameters  for  the  various  gases  and  sources.  It  is  readily  observed  that  xg  approaches  xs  as 
P0  increases  or  T0  decreases,  and  supercooling  is  seen  to  range  approximately  from  30  to 
200  K. 

4.2  DESCRIPTION  OF  SCALING  ANALYSIS 

Empirical  functional  relations  of  f  with  x,  xg,  P o)  T0,  and  D  or  D,  were  obtained  (Refs.  22 
and  23)  using  the  results  of  the  axial  variation  of  the  scattering  function  data.  These 
variations  in  the  condensate  growth  region  were  represented  by 


and  by  graphical  determination  it  was  found  that 

b  cc  P;m°  T“t  (8) 

se  -  (p"!  D«,  T”‘  <9> 

The  quantity  Deq  is  D  and  (D,  cot  fl'/JQ?)  for  sonic  orifices  and  conical  nozzles, 
respectively.  The  expression  C(y)  is  a  parameter  dependent  upon  specific  heat  ratio  and  is 
defined  in  Ref.  24.  The  sonic  orifice  scaling  constants  are  tabulated  in  Table  2,  and  the 
conical  nozzle  constants  are  given  in  Table  3. 

By  using  the  well-depth  and  range  parameters,  e  and  a,  respectively,  of  the  Lennard- 
Jones  12-6  intermolecular  potential  function,  reduced  onset  pressures  (Pg)  and 
temperatures  (Tg)  have  been  determined  using 

P*e  =-  P$/U  a3)  (10) 
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t;  =  T0/(f/k)  (id 

where  k  is  Boltzmann’s  constant.  The  reduced  onset  parameters  are  also  given  in  Table  1, 
and  the  intermolecular  potential  constants  which  were  used  are  given  in  Table  4.  The  loci  of 
condensation  onset  are  shown  in  Figs.  34  through  37.  The  common  locus  for  the 
homonuclear  diatomic  molecules  N2  and  O2  is  observed.  In  addition,  the  difference  between 
the  Ar  results  and  the  results  for  N2  and  O2  is  noted  to  illustrate  the  effect  of  specific  heat 
ratio  upon  the  location  of  the  onset  locus.  Not  unexpectedly,  the  loci  of  polar  diatomics  such 
as  CO  and  HC1  are  found  to  be  substantially  different,  and  this  illustrates  the  inadequacy  of 
the  two-parameter  Lennard-Jones  potential  for  describing  the  interaction  of  polar 
molecules.  As  is  well  known,  the  interaction  of  polar  molecules  includes  important 
contributions  from  dipole-dipole  forces,  and  the  three-parameter  Stockmayer  potential  is 
more  accurate  for  this  description. 

4.3  CORRELATION  OF  STATIC  PROPERTY  MEASUREMENTS,  RAYLEIGH 

SCATTERING  MEASUREMENTS,  AND  THEORETICAL  PREDICTIONS 

WITH  CONDENSATION 

Figures  38  and  39  show  the  measured  axial  variation  of  I'(!),  n(N2)/nD,  and  T/Tc  for 
two  P0  values  for  N2  nozzle  flows.  Number  density  and  temperature  were  obtained  using 
laser-Raman  scattering  and  also  by  electron  beam  fluorescence.  A  more  complete  exposition 
of  the  static  parameter  measurements  of  N2  and  other  gas  species’  flows  can  be  found  in 
Refs.  5,  14,  18,  and  19,  but  the  measurements  shown  in  Figs.  38  and  39  well  illustrate  the 
results  of  the  measurements.  As  shown  in  Figs.  38  and  39,  the  monomer  number  density  is 
little  affected  by  the  condensation  process;  however,  the  static  temperature  can  be  increased 
by  as  much  as  50  per  cent  above  the  isentropic  prediction  due  to  the  heat  release  in  the 
condensation  process.  It  is  also  observed  that  the  increase  in  temperature  correlates  very  well 
with  the  onset  of  condensate  growth.  Not  shown  in  either  Fig.  38  or  39  are  the  far-field 
measurements  of  n(N2)/n0  made  with  the  electron  beam  system;  however,  these  results  are 
given  in  Ref.  4  and  show  excellent  agreement  with  the  MOCS  prediction. 

Predictions  of  Rayleigh  scattering  intensity,  condensate  mass  fraction,  and  static 
temperature  generated  using  the  liquid  drop,  monodisperse  condensation  model  are  also 
shown  in  Figs.  38  and  39,  as  well  as  in  Fig.  19.  This  model  is  described  in  detail  in  Ref.  18. 
Briefly,  this  calculation  assumes  the  condensing  flow  field  to  be  inviscid,  adiabatic,  and  one¬ 
dimensional,  with  no  mass  transfer  across  the  stream  tube  boundary.  The  gas  and 
condensate  are  assumed  to  obey  the  perfect  gas  relation.  The  condensed  phase  is  assumed  to 
be  of  the  form  of  monodisperse  spherical  drops  or  particles  which  are  characterized  by  bulk 
properties  and  to  be  in  the  free  molecular  flow  regime  relative  to  the  uncondensed  phase. 
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Condensate-gas  velocity  slip  effects  are  ignored,  and  the  condensate  growth  rate  is 
determined  by  gas-condensate  interaction  only.  The  mass  accommodation  coefficient  is 
assumed  to  be  unity.  Initial  size  and  number  density  of  spontaneous  nucleation  sites  are 
adjustable  parameters  and  may  be  selected  so  as  to  reproduce  as  closely  as  possible  the 
experimental  Rayleigh  scattering  results.  For  the  predictions  shown  in  Figs.  38,  39,  and  19 
the  initial  nucleation  sites  were  chosen  to  be  dimers  with  mole  fractions  on  the  order  of  10*3, 
and  the  starting  point  for  the  calculation  was  the  saturation  point,  xs.  The  model  can  make 
adequate  (±50  percent)  predictions  for  pure  gas  expansions  of  the  Rayleigh  scattering 
intensity  variation  and  the  location  of  the  condensation  onset  if  the  initial  nucleation  sites 
are  assumed  to  be  dimers  with  a  concentration  equal  to  the  equilibrium  dimer  mole  fraction 
at  saturation.  The  predictions  of  the  model  are  observed  to  agree  more  closely  with  the 
experimental  measurements  as  degree  of  condensation  increases.  As  can  be  best  observed  in 
Fig.  38,  the  model-predicted  temperature  increase  that  results  from  the  condensation  process 
agrees  well  with  the  experimentally  determined  temperatures.  Condensate  mass  fractions  on 
the  order  of  10~2  to  10-'  are  predicted  for  the  cases  shown  in  Figs.  38,  39,  and  19.  As  noted  in 
Figs.  38  and  39  the  number  of  molecules  per  cluster  is  100,  and  the  size  of  the  clusters  is  on 
the  order  of  10  A. 

Also  shown  in  Figs.  25  and  38  are  measurements  of  the  Rayleigh  scattering 
depolarization  ratio,  p,  which  is  defined  as 

p  =  ru)/p(!)  (12) 


The  axial  value  of  p  (Fig.  38)  does  indeed  decrease  rapidly  from  its  room  temperature 
monomer  value  as  the  cluster  growth  region  is  entered.  Furthermore,  the  radial  value  of  p 
(Fig.  25)  decreases  rapidly  as  the  off-axis  cluster  growth  region  is  entered.  This  behavior  is 
intuitively  expected,  because  as  the  linear  N2  molecules  cluster  together  they  should  form 
more  spherical  scatterers  that  contribute  a  larger  portion  to  the  scattered  intensity  than  do 
the  monomers. 

5.0  DISUCSSION  OF  THRUSTER  AND  SIMULATED  THRUSTER  DATA 
5.1  THRUSTER  RESULTS 

The  large  quantity  of  data  accumulated  during  the  thruster  experiments  is  presented  in 
detail  in  Refs.  16  and  25;  therefore,  only  data  pertinent  to  the  phenomenon  of  condensation 
will  be  discussed  here.  Rayleigh  scattering  levels,  I' (II),  measured  in  the  aged  and  new 
thruster  plumes  are  compared  in  Table  5  for  two  test  conditions  (2A  and  2S;  see  Ref.  16). 
The  scattering  function,  f,  was  evaluated  for  a  7  =  1.2  plume  expansion  and  is  tabulated  in 
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Table  5.  As  observed  in  Table  5,  the  f  factors  for  the  aged  thruster  are  approximately  a 
factor  of  50  greater  than  those  for  the  new  thruster.  It  is  believed  that  the  high  Rayleigh 
scattering  levels  for  the  aged  thruster  are  a  result  of  raw  hydrazine  traveling  down  the  axial 
centerline.  This  belief  is  supported  by  the  high  mass  deposition  rate  measured  on  the  axial 
centerline  by  a  QCM,  particle  sampling  measurements,  and  radial  Rayleigh  scattering 
profiles  which  reveal  the  extremely  high  levels  only  on  the  axial  centerline  (see  Ref.  16). 

The  Rayleigh  scattering  levels  for  the  new  thruster  are  an  order  of  magnitude  or  more 
greater  than  expected  for  a  noncondensing,  particulate-free  plume  expansion  as  can  be 
observed  in  Tables  5,  6,  and  7  and  in  Fig.  40,  which  demonstrates  the  axial  variation  of  I '  (II) 
for  the  new  thruster.  It  can  be  seen  in  Fig.  40  that  the  measured  values  of  I  '(II)  exceed  the 
predicted  values  by  a  magnitude  which  increases  with  increasing  x,  and  this  is  well  illustrated 
by  the  axial  variation  of  the  scattering  function,  f,  plotted  in  Fig.  41  for  thruster  test 
conditions  2S  and  2 C.  Both  of  these  conditions  possess  an  axial  coordinate  dependence 
which  exhibits  a  rapid  increase  at  the  lower  x  values  and  becomes  asymptotic  at  f  =  100  for 
the  far-field  plume.  Although  the  l'(li)  values  for  the  thruster  were  obtained  using 
calibration  factors  appropriate  for  pure  N2,  the  analysis  in  Ref.  25  has  shown  that  f  should 
be  zero  with  an  uncertainty  on  the  order  of  ±  0.1  for  a  noncondensing  plume  expansion  of 
N2,  H2,  and  NH3  with  1.2  <  7  <  1.3.  The  possible  presence  of  either  particulate  matter  or 
raw  fuel  droplets  cannot  explain  the  Rayleigh  scattering  results,  for  if  such  flow-borne 
particulates  followed  the  expansion,  then  f  would  be  approximately  constant;  if  the  material 
were  concentrated  along  the  axial  streamtube,  then  f  would  increase  continually  with  x. 

The  Raman  measurements  of  temperature  and  total  density  in  the  plume  are  also  shown 
in  Fig.  40;  assuming  the  expansion  to  have  frozen  7  and  chemistry  and  using  the  isentropic 
relation  between  temperature  and  total  number  density,  a  plume  7  of  1.22  is  determined 
(Ref.  25).  The  Raman  and  mass  spectrometric  measurements  of  species  mole  fractions 
enable  the  inferential  determination  of  a  plume  7  of  approximately  1 .30  (Ref.  25).  These 
results  demonstrate  that  the  7  =  1.3  expansion  curve  should  be  followed  but  that  7  =  1.22  is 
the  actual  heat  capacity  ratio.  This  situation  is  indicative  of  significant  condensation  in  the 
flow  field,  and  Fig.  13  of  Ref.  18  shows  a  similar  expansion  history  and  behavior  regarding 
N2  condensation.  Therefore,  the  Raman  and  mass  spectrometric  measurements  as  well  as  the 
measured  magnitude  and  variation  of  the  Rayleigh  scattering  function  support  the 
contention  that  significant  condensation  is  occuring  in  the  thruster  plume. 

Using  the  results  tabulated  in  Tables  6  and  7,  one  can  scale  I' (11)  as  a  function  of 
reservoir  pressure  and  temperature  for  an  initial  catalyst  bed  temperature  of  478  K  and  for  x 
=  45.2.  The  scaling  is 

I'(l)  «p*-5t;5 
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with  a  correlation  coefficient  of  -0.994.  Unfortunately,  a  change  in  initial  catalyst  bed 
temperature  with  the  average  reservoir  parameters  remaining  the  same  destroys  the  scaling 
as  shown  in  Fig.  42.  Therefore,  it  is  not  possible  to  use  only  the  reservoir  parameters  P0  and 
T0  for  scaling  of  I '(II)  in  the  thruster  plume. 

5.2  SIMULATED  THRUSTER  RESULTS 

In  an  additional  effort  to  investigate  the  large  levels  of  Rayleigh  scattering  in  the  thruster 
plume,  the  thruster  plume  expansion  was  simulated  in  the  RVC.  A  nozzle  assembly  was  used 
that  duplicated  the  nozzle  of  the  thruster,  and  a  set  of  three  gas  mixtures  was  prepared  by  the 
AEDC  Chemical  Laboratory  to  simulate  the  thruster  plume  composition.  The  first  two 
mixtures  were  prepared  to  simulate  a  plume  y  of  1.2,  and  the  third  mixture  simulated  a 
plume  7  of  1.25.  The  first  mixture  was  a  binary  mixture  of  Nz  and  NH3.  The  NH3  mole 
fraction  was  equal  to  that  predicted  by  CONTAM  II  for  the  thruster  test  condition  2S,  and 
the  Ni  mole  fraction  was  equal  to  the  sum  of  the  H2  and  Ni  mole  fractions  predicted  by 
CONTAM  II  for  the  2S  condition.  For  the  second  mixture,  the  N2,  H2,  and  NH3  mole 
fractions  were  those  predicted  by  CONTAM  II  for  test  condition  2S.  The  mole  fractions  of 
N2,  H2,  and  NH3  for  the  third  mixture  were  determined  using  Fig.  90  of  Ref.  16. 

The  mixtures  were  permitted  to  flow  steadily  through  a  heated  source/nozzle  assembly, 
and  the  reservoir  pressure  and  temperature  were  set  to  equal  the  values  of  the  2S  condition  as 
given  in  Ref.  25.  The  flow  could  be  maintained  for  approximately  60  sec  when  H2  was  a 
mixture  constituent  before  the  RVC  pressure  reached  a  level  of  10  mtorr. 

The  results  of  these  experiments  are  shown  in  Fig.  40  and  are  listed  in  Table  5  for 
x  =  78.5.  It  can  be  observed  that  the  7  =  1.2  and  1.25  mixutres  of  N2-H2-NH3  have 
Rayleigh  scattering  levels  that  bracket  those  observed  for  the  new  thruster.  The  f  factor  for 
the  7  =  1.2,  N2-H2-NH3  mixture  is  approximately  a  factor  of  2  higher  and  the  f  factor  for 
the  7  =  1.25,  N2-H2-NH3  mixture  is  approximately  a  factor  of  2  lower  than  the  f  factor  for 
the  new  thruster  at  the  2S  condition.  An  axial  variation  of  Rayleigh  scattering  for  the  7  = 

1 .2,  N2-H2-NH3  mixture  is  also  shown  in  Fig.  40,  and  it  can  be  observed  that  the  behavior  is 
very  similar  to  that  observed  for  condensation  in  the  pure  gas  and  binary  mixture  expansions 
in  the  scaling  law  experiments.  As  x  decreases,  the  Rayleigh  scattering  is  observed  to 
approach  the  level  expected  for  an  isentropic,  noncondensing  expansion.  It  is  speculated 
that  the  same  behavior  would  have  been  observed  for  the  thruster  if  reliable  measurements 
could  have  been  made  for  x  <  20. 

Analysis  of  the  intersection  of  N2H4,  H2O,  NH3,  and  N2  vapor  pressure  curves  with  7  = 

1.2  and  1.3  isentropes  for  test  condition  2S  reservoir  conditions  (Ref.  25)  shows  that  NH3 
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saturation  will  occur  from  x  =  9.5  to  34  and  that  H20  and  N2H4  will  saturate  much  earlier. 
The  previous  experiments  with  N2/C02  mixtures  indicate  that  the  H20  and  N2H4 
condensate  will  provide  nucleation  sites  that  will  bring  about  NH3  condensation  near  the 
saturation  point.  This  would  explain  the  lack  of  a  pronounced  peak  in  the  Rayleigh 
scattering  axial  profiles  for  the  thruster  and  the  presence  of  the  peak  for  the  simulated 
mixture  profile.  It  is  therefore  believed  that  the  high  levels  of  Rayleigh  scattering  in  the  new 
thruster  plume  are  a  direct  result  of  condensation  of  gas  mixture  species  in  the  flow  field. 

6.0  SUMMARY 


6.1  CONCLUSIONS 


A  multi-year  experimental  program  for  characterization  of  the  onset  and  growth  of 
condensation  in  expansion  flow  fields  with  regard  to  nozzle  geometry,  reservoir  conditions, 
molecular  parameters,  and  flow  composition  has  been  completed.  This  report  has  reviewed 
the  experimental  results;  a  number  of  conclusions  can  be  made  and  are  listed  as  follows: 

1 .  The  scaling  law  for  the  pure  gas  scattering  function  can  be  written  as 

f «  T^mi  Df2/cot  6Yi  (Conical  Nozzles) 

f  «  P™°  Dm2  (Sonic  Orifice) 


2.  The  scaling  law  for  axial  centerline  location  of  condensation  onset  for  pure 
gases  can  be  written  as 


oe 


Tm4 


3.  Intermolecular  parameters  were  used  for  normalization  of  the  pure  gas 
condensation  onset  pressure  and  temperature  values,  and  the  condensation 
onset  loci  were  plotted  in  the  P*-T*  plane.  These  loci  represent  the  practical 
limit  to  which  a  pure  gas  can  be  expanded  (supercooled)  before  massive 
homogeneous  condensation  begins.  The  loci  are  as  follows: 

Pjg  -  7.12  x  10'3  Tg4  24  Sonic  Orifice  \ 

f  Argon 

P^  =  2.10  x  10'3  Tg3,66  Conical  Nozzle  ' 
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V*q  -  3.68  >  10'2  Tg5'9'1 
Pq  =  1.02  *  IO'aT5s-16 
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4.  The  condensation  onset  loci  results  imply  that  classes  of  molecules  well 
described  by  the  12-6  Lennard- Jones  potential  will  have  common  onset  loci 
depending  upon  the  specific  heat  ratio.  However,  polar  molecules  better 
represented  by  the  12-6-3  Stockmayer  potential  do  not  show  common  onset  loci 
for  common  specific  heat  ratios. 


5.  A  liquid-drop,  monodisperse  distribution  condensation  calculation  for  flow 
fields  of  pure  gases  was  developed.  This  model  can  make  ±  50-percent 
predictions  of  the  axial  variation  of  observed  Rayleigh  scattering  by  choosing 
the  initial  nucleation  sites  to  be  dimers  and  by  choosing  the  initial  dimer 
concentration  to  be  the  equilibrium  dimer  concentration  at  the  saturation 
point.  Furthermore,  the  model  predicts  the  observed  static  temperature  rise 
caused  by  condensation  within  ±20  percent.  Cluster  size  and  condensate  mass 
fraction  are  also  predicted  by  the  model. 


6.  Rayleigh  scattering  radial  profiles  of  pure  gas  flows  have  shown  that  off-axis 
condensation  has  significantly  preceded  the  axial  phenomenon. 
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7.  Rayleigh  scattering  axial  profiles  of  binary  gas  mixtures  (N2/CO2  and  NH3/N2) 
indicate  that  the  nucleation  of  the  easily  condensable  species  (CO2  and  NH3) 
will  provide  nucleation  sites  that  will  cause  N2  condensation  near  the  saturation 
point  rather  than  permitting  the  large  degrees  of  supercooling  observed  for  the 
pure  gas  expansion. 

8.  Rayleigh  scattering  measurements  as  well  as  Raman  scattering  and  mass 
spectrometric  measurements  indicate  a  significant  quantity  of  condensate  in  the 
new  (refurbished)  thruster  plume,  and  this  was  subsequently  verified  by 
Rayleigh  scattering  measurements  in  a  particulate-free  simulated  thruster 
plume. 

9.  The  Rayleigh  scattering  intensity  from  the  new  thruster  plume  scales  as  Pq2-5T '5 
for  a  given  initial  catalyst  bed  temperature  in  dramatic  contrast  to  the  scaling 
observed  for  pure  gases.  Because  of  a  significant  dependence  upon  initial 
catalyst  temperature,  Rayleigh  scattering  intensity  cannot  be  scaled  using  only 
reservoir  parameters  P0  and  T0. 

10.  Monopropellant  thruster  plumes  can  be  successfully  simulated  using  gas 
mixtures  for  the  purpose  of  studying  condensation  phenomena  in  the  plume. 

6.2  RECOMMENDATIONS  FOR  FUTURE  WORK 

The  primary  objective  for  future  work  concerning  flow-field  condensation  should  be  to 
develop  an  accurate  computer  code  for  predicting  the  effects  of  condensation  in  expansion 
flows.  It  is  recommended  that  the  development  of  the  code  proceed  in  the  following  fashion. 

a.  The  present  liquid-drop,  monodisperse  condensate  growth  model  developed  at 
AEDC  for  predictions  of  pure  gas  axial  centerline  condensation  effects  should 
be  expanded  for  the  purpose  of  prediction  of  off-axis  condensation.  The  results 
should  be  compared  to  the  experimental  results  of  this  report. 

b.  The  classical  capillarity  theory  for  calculation  of  the  number  density  of  initial 
nucleation  sites  and  their  size  should  be  added  to  the  condensate  growth  model 
developed  at  AEDC.  The  results  should  be  compared  to  the  experimental  results 
of  this  report  to  quantitatively  assess  the  shortcomings  of  the  classical 
condensation  theory. 
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c.  An  alternate  approach  for  calculation  of  the  number  density  and  size  of  initial 
nucleation  sites  should  be  developed,  and  it  should  be  strictly  kinetic  in  nature. 
That  is,  the  kinetic  formation  of  dimers,  trimers,  etc.  must  be  calculated  and 
subsequently  used  with  the  condensate  growth  model.  Again,  the  results  should 
be  compared  to  the  experimental  results  of  this  report  to  assess  the  validity  of  the 
model. 
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Figure  1 .  Experimental  arrangement  for  the  scaling  law 
and  simulated  thruster  studies. 
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Figure  2.  Sonic  orifice  and  conical  nozzle  schematics. 
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Figure  3.  Thruster  installation  in  the  RVC. 


Figure  4.  Experimental  arrangement  for  special  diagnostic 
instrumentation  for  thruster  studies. 
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Figure  5.  Axial  variation  of  l'(!l)  for  all  N2 
reservoir  pressures  and  sonic 
orifice  diameters  investigated. 


Figure  6.  Axial  variation  of  l'(ll)  for  all  Oz 
reservoir  pressures  investigated, 
sonic  orifice. 
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Figure  7.  Axial  variation  of  l'(ll)  for  all  CO 
reservoir  pressures  investigated, 
sonic  orifice. 
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Figure  11.  Axial  variation  of  l'( II)  for  all  N2 
reservoir  pressures  investigated, 
0>s  =  14.5  deg. 


Figure  12.  Axial  variation  of  I' ( II )  for  all  02 
reservoir  pressures  investigated, 
054  =  10.5  deg. 
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Figure  13.  Axial  variation  of  l'(ll)  for  all  CO 
reservoir  pressures  investigated, 
0%  =  10.5  dea. 
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Figure  14.  Axial  variation  of  !'( H)  for  ail  NO 
reservoir  pressures  investigated, 
0%  =  10.5  deg. 
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Figure  16.  Axial  variation  of  l'(ll)  for  all  H20  reservoir 
pressures  investigated,  0u  =10.5  deg. 
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Figure  17.  Axial  variation  of  I' (II)  for  all  C02 
reservoir  pressures  investigated, 

9ji  =  10.5  deg. 


Figure  18.  Axial  variation  of  l'(ll)  for  all  Ar 
reservoir  pressures  investigated, 
0j4  =  10.5  deg. 
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Figure  20.  Axial  variation  of  l'(ll)  for  all  02 
reservoir  temperatures  investigated, 
6%  -  10.5  deg. 
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Figure  21.  Axial  variation  of  I' (II)  for  all  CO  Figure  22.  Axial  variation  of  l'( II)  for  all  HCI 

reservoir  temperatures  investigated,  reservoir  temperatures  investigated, 

0j4  =  10.5  deg.  054  =  10.5  deg. 
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Figure  26.  Radial  variation  of  l'(ll)  for  an  N2 
reservoir  pressure  of  7.85  atm  at 
six  axial  positions,  T0  =  288  K, 
0j4  =  10.5  deg. 


Figure  27.  Radial  variation  of  i'(!l)  for  an  N2 
reservoir  pressure  of  7.85  atm  at 
six  axial  positions,  T0  =  355  K, 
6%  =  105  deg. 
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Figure  32.  Axial  variation  of  scattering 

function,  f,  for  10.5-deg  nozzle, 
N2  gas,  T0  variation. 
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Figure  33.  Diagram  of  expansion  process. 
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Figure  34.  N2  and  02  condensation  onset  locus  using 
reduced  values. 


Figure  35.  Ar  and  CO  condensation  onset  locus  using 
reduced  values. 


AEDC-TR-80-16 


o  nio 

A  ttyNo 
*  WTo 

a  Depolarization  Ratio,  p 


P0*  10. 2  atm,  N2 

14. 5 -deg  Half-Angle  Conical  Nozzle 


10°  101  102  103 


Figure  38.  Axial  variation  of  l'(ll)  number  density, 
temperature,  mass  fraction  (g),  and 
depolarization  ratio:  14.5-deg  conical 
nozzle,  P0  =  10.2  atm,  N2 . 
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Figure  39.  Axial  variation  of  l'(ll),  number  density, 
temperature,  and  mass  fraction:  14.5-deg 
conical  nozzle,  P0  =  6.80  atm,  N2. 
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41.  Axial  variation  of  the  scattering  function,  f, 
in  the  new  thruster  plume. 
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Figure  42.  Variation  of  Rayleigh  scattering  in  the  new  thruster  plume 
with  reservoir  parameters  at  axial  position  x  =  45.2. 
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Table  1 


Ol 


O 


ndensation  Onset  Parameters 


V 

torr 

V 

K 

p*e 

T$ 

<v 

«*>*■ 

K 

10.8 

0.29 

16.9 

1.49 

X 

10-6 

0.178 

71.1 

41.1 

12. 5 

0.125 

14.9 

6.44 

X 

10'7 

0.157 

103.3 

40.6 

14.2 

0.063 

13.3 

3.25 

X 

io"7 

0.140 

125.4 

39.7 

16.4 

0.0314 

11.9 

1.62 

X 

ID'7 

0.L25 

156.0 

39.1 

9.64 

0.274 

18.2 

1.41 

X 

ID’6 

0.191 

53.3 

38.2 

11.75 

0.110 

15.1 

5.67 

X 

IQ'7 

0.159 

82.7 

38.9 

21.9 

0.162 

15.1 

8.35 

X 

to'7 

0.159 

101.9 

40.9 

27.7 

0.045 

11.7 

2.32 

X 

io-7 

0.123 

209.0 

42.3 

31.5 

0.0176 

9.9 

9.07 

X 

io-8 

0.104 

329.5 

42.1 

39.5 

0.0064 

8.0 

3.30 

X 

io'8 

0.084 

563.0 

42.0 

22.6 

0.459 

18.8 

2.37 

X 

io-6 

0.198 

49.1 

38.7 

24.8 

0.275 

17.4 

1.42 

X 

10"6 

0.183 

56.7 

38.6 

27.2 

0.16 

15.6 

8.25 

X 

io-7 

0.164 

31.3 

39.4 

31.4 

0.0579 

12.9 

2.98 

X 

10-7 

0.136 

140.0 

40.1 

37.0 

0.0222 

10.6 

1.14 

X 

10-n 

0.112 

250.0 

40.9 

41.0 

0.0106 

9.5 

5.46 

X 

io-8 

0.0999 

323.0 

40.0 

24.4 

0.309 

17.02 

1.59 

X 

io-6 

0.179 

67.3 

39.7 

32.4 

0.0990 

13.89 

5.10 

X 

io'7 

0.146 

121.2 

40.8 

40.0 

0.0461 

12.08 

2.38 

X 

10'7 

0.127 

180.9 

41.1 

50.0 

0.0211 

11.04 

1.10 

X 

io-7 

0.116 

218.0 

40.3 

22.4 

1.19 

23.9 

6.13 

X 

io-6 

0.251 

22.7 

34.6 

27.6 

0.40 

19.1 

2.06 

X 

io-6 

0.201 

45.5 

37.5 

31.5 

0.138 

15.4 

7.11 

X 

io-7 

0.162 

81 .9 

39.4 

37.8 

0.046 

12.5 

2.37 

X 

io-7 

0.132 

152.2 

40.3 

42.0 

0.0235 

11.2 

1.21 

X 

io-7 

0.118 

206.0 

40.1 

— 

— 

— 

— 

— 

— 

— 
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Table  1.  Continued 
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Table  1.  Continued 


to 


pe 

TS 

(S0>” 

K 

2.68  x  Uf6 

23.2 

45.8 

1 .19  x  10"6 

30.9 

46.5 

5.61  x  10"7 

0.114 

43.9 

46.4 

1.78  x  10-7 

0.0798 

96.5 

49.4 

2.6  x  10"6 

0.155 

21.2 

44.9 

1.32  x  10-6 

0.131 

29.4 

46.6 

8.6  x  10"7 

0.119 

34.7 

46.9 

2.09  x  10~7 

0.0778 

96.6 

50.6 

1.06  x  10-6 

0.124 

32.4 

47.2 

5.38  x  10"7 

0.110 

42.2 

46.9 

m 

i 

o 

X 

o 

m 

0.338 

23.66 

103.98 

9.43  x  10-6 

0.313 

30.12 

107.03 

5.68  x  10"6 

0.29 

38.36 

108.67 

3.21  x  10-6 

0.264 

48.35 

111.27 

1.35  x  10“6 

0.231 

69.67 

113.87 

3.81  x  10-6 

0.266 

50.14 

113.47 

2.49  x  10"6 

0.258 

52.23 

111.1 

1 .09  x  10"7 

0.149 

27.3 

78.2 

2.35  x  10'7 

0.150 

33.3 

84.7 

7.26  x  10-7 

0.164 

29.7 

88.2 

1.03  x  10-6 

0.165 

31 .0 

90.9 

1.92  x  10"6 

0.179 

25.6 

90.2 

6.02  x  10-7 

0.163 

35.8 

88.7 

5.86  x  10"7 

0.167 

25.7 

84.1 

4.41  x  10"7 

0.178 

17.7 

75.7 
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Table  1.  Concluded 


P..  *• 


(s.r  (Be)* 

8  K 


4.25 

4.7 

92 

9.05 

0.805 

54.9 

1.89  x  10'6 

0.419 

5.84 

37.1 

4.0 

10.7 

95 

4.89 

6.9 

83.8 

1.62  x  10-5 

0.64 

1.55 

11.2 

3.9 

17.8 

97.5 

— 

— 

— 

— 

— 

— 

— 

3.8 

24.6 

99 

— 

— 

— 

— 

— 

— 

— 

3.75 

31.0 

100 

— 

— 

— 

— 

— 

— 

— 

3.7 

39.5 

101 

— 

— 

— 

— 

— 

— 

— 

0.178 

16.9 

290 

12.2 

0.0902 

81.5 

4.28  x  10"8 

0.215 

187.4 

208.5 

43.5  309 

61.5  315 

41.0  308 

21.0  295 

360  354 

31.0  301 

29.5  300 

28.2  300 

25.5  298 

23.5  296 


78.0  320 
45.1  309 
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Table  2.  Sonic  Orifice  Scaling  Constants 


Gas 

D, 

mm 

m 

0 

m2 

m3 

N2 

1.325 

2.94 

1.84 

0.25 

N2 

3.050 

2.93 

1 .84 

0.25 

°2 

1.325 

2.98 

2.05 

— 

CO 

1.325 

2.87 

— 

— 

Ar 

3.200 

3.03 

1 .96 

— 

Table  3.  Conical  Nozzle  Scaling  Constants 
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Table  4.  Lennard-Jones  12-6  Intermolecular  Potential  Constants 


Gas 

T'  =  e/k, 

K 

ergs 

D'  =  a, 
cm 

P"  =  e/cf3, 
torr 

N2 

95.05 

1.31  x  10"U 

3.70  x  10"8 

1.94  x  105 

°2 

117.50 

1.62  x  10"14 

3.58  x  10~8 

2.65  x  105 

CO 

100.20 

1.38  x  10“U 

3.76  x  10-8 

1.94  x  105 

C°2 

205.00 

2.83  x  10"U 

4.07  x  10-8 

3.15  x  105 

Ar 

119.80 

1 .65  x  10"14 

3.40  x  10"8 

3.14  x  105 

NO 

131.00 

1.81  x  10"14 

3.17  x  10"8 

4.26  x  105 

HC1 

328.00 

-14 

4.53  x  10 

3.36  x  10"8 

8.96  x  105 

h2o 

380.00 

5.24  x  10"14 

2.65  x  10-8 

2.11  x  106 
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Table  5.  Comparison  of  Rayleigh  Scattering  Data  for  the  Aged, 
New,  and  Simulated  Thruster  Plumes 


Plume 

Test 

Condition 

I'd  1) 

MOCS 

n/nQ,  Y  =  1-2 

MOCS 

n/n  ,  y  =  1 *3 
o 

Aged 

Thruster 

2S 

1.04  x  10"1 

5.40  x  10“5 

1.1  x  10_A 

New 

Thruster 

2S 

2.08  x  10“3 

Simulated 

Thruster, 

Y  =  1.2 

2S, 

n2-nh3 

6.06  x  10“3 

Simulated 

Thruster, 

2S 

n2-nh3:h2 

3.94  x  10"3 

Y  =  1.2 

Simulated 

2S 

n2-nh3-h2 

Thruster, 

Y  -  1.25 

1.29  x  10 

Aged 

Thruster 

2A 

4.06  x  10_1 

1.03 

K  lO"3 

1.3  x  10-3 

New 

Thruster 

2A 

8.14  x  10“3 

1.03 

K  10“3 

1.3  x  10"3 

f  (Y  =  1.2) 

1.92  x  103 
3.75  x  101 

1.11  x  102 
7.19  x  101 

2.29  x  101 

3.93  x  102 


New 

Thruster 


2A 


6.90 
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Table  6.  Rayleigh  Scattering  Intensity,  New  Thruster  Plume, 
Test  Period  No.  25,  x  =  45.2 


Test 

Condition 

I'(ll) 

To,K 

-3 

n  jCm 
o 

Y  = 
n/n 

1.2 

o 

2S 

3.71  x  10“3 

653 

7.00  x  1019 

4.3  x 

10-4 

2A 

3.61  x  10“3 

710 

9.67  x  1019 

2B 

3.04  x  10“3 

739 

1.186  x  1020 

2C 

2.78  x  10“3 

757 

1.326  x  1020 

1 A 

5.07  x  10"3 

679 

20 

1.013  x  10 

1C 

2.49  x  10-3 

730 

20 

1.372  x  10  u 
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Table  7.  Rayleigh  Scattering  Intensity,  New  Trhuster  Plume, 
Test  Period  No.  26 


A 

X 

Test 

Condition 

2S 

Test 

Condition 

2A 

Test 

Condition 

2B 

Test 

Condition 

2C 

Test 

Condition 

1A 

Test 

Condition 

1C 

278.5 

3.35  x  10~4 

3.55  x  10"4 

2.32  x  10“4 

1.96  x  10~4 

2.48  x  10~4 

2.17  x  10"4 

211.8 

5.34  x  10"4 

4.12  x  10“4 

3.71  x  10"4 

3.80  x  10~4 

4.11  x  10"4 

2.35  x  10~4 

145.2 

8.29  x  10~4 

6.38  x  10~4 

5.97  x  10~4 

4.81  x  10~4 

6.51  x  10"4 

3.95  x  10“4 

111.8 

1.57  x  10"3 

8.47  x  10“4 

1.26  x  10"3 

9.67  x  10-4 

1 .28  x  10-3 

5.95  x  10"4 

78.5 

2.08  x  10~3 

1.21  x  10”3 

1.50  x  10-3 

1.68  x  10~3 

1.99  x  10“3 

1.17  x  10-3 

-3 

-3 

45.2 

3.83  x  10 

3.10  x  10 

— 

— 

— 

— 

28.5 

1 .12  x  10-2 

8.14  x  10"3 

8.15  x  10-3 

5.75  x  10-3 

7.98  x  10"3 

4.80  x  10"3 
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NOMENCLATURE 

a  Radius  of  scatterer 

C  (  7  )  7  -dependent  parameter 

CONTAM  II  Computer  code  for  predicting  plume  contamination  from  liquid 
monopropellant  and  bipropellant  rocket  engines  on  spacecraft  surfaces 

D  Sonic  orifice  diameter 

Dt  Conical  nozzle  throat  diameter 

Deq  Equivalent  diameter  defined  as  D  for  sonic  orifices  and  Cfy)  Dt  cot 0Vl 

for  conical  nozzles 

f  Rayleigh  scattering  function  defined  by  Eq.  (4) 

g  Condensate  mass  fraction 

I  Relative  Rayleigh  scattering  intensity  defined  by  Eq.  (1) 

I0  Laser  beam  intensity 

I' (II),  I'(-L)  Relative  Rayleigh  scattering  intensity,  normalized  to  the  relative 
Rayleigh  scattering  intensity  of  a  gas  sample  of  number  density  n0, 
polarized  parallel  and  perpendicular  to  the  plane  of  polarization  of  the 
incident  laser  beam,  respectively 

Constant  in  Eq.  (1) 

Boltzmann’s  constant 

Method  of  characteristics  solution 

Scaling  constants 

Number  density  of  gas  species  and  i-mer  number  density,  respectively 
N2  number  density 
Reservoir  number  density 
Local  total  number  density 
Reservoir  pressure 


K 

k 

MOCS 

m0,mi,m2, 

1113,014 

n,  nj 

n(N2) 

n0 

nt 

Po 
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Ps 


P  v(? 

Pg,  p; 


QCM 

RVC 

r 

(s)° 

(S  g)° 

T 

To 

Ts 

t„,t; 


x;x 


xs 


a,aj 


7 

€ 

0i/2 

Xo 

P 

a 


Saturation  pressure 

Equilibrium  vapor  pressure  at  condensation  onset 

Pressure  at  condensation  onset  and  Pg  normalized  by  the  intermolecular 
potential  parameter  e/ <r3 

Quartz  crystal  microbalance 

Research  Vacuum  Chamber 

Radial  distance  from  flow-field  centerline 

Isentropic  supersaturation 

Isentropic  degrees  of  supercooling 

Static  temperature 

Reservoir  temperature 

Saturation  temperature 

Temperature  at  condensation  onset  and  Te  normalized  by  the 
intermolecular  potential  parameter  e/k 

Axial  position  in  the  flow  field;  normalized  axial  position  in  the  flow 
field  (x  =  x/D  for  sonic  orifices  and  x  =  x/D,  for  conical  nozzles) 

Axial  location  of  saturation 

Axial  location  of  condensation  onset 

Electronic  polarizability  and  i-mer  electronic  polarizability,  respectively 
Specific  heat  ratio 

Well  depth  of  Lennard-Jones  intermolecular  potential  function 
Expansion  half-angle  of  conical  nozzle 
Wavelength  of  incident  laser  beam 
Depolarization  ratio  of  Rayleigh  scattering 

Range  parameter  of  Lennard-Jones  intermolecular  potential  function 
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